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Multidrug and Toxic compound Extrusion proteins (MATE) are a group of secondary active transporters 
with ubiquitous occurrences in all domains of life. This is a newly characterized transporter family with 
limited functional knowledge in plants. In this study, we functionally characterised two members of rice 
MATE gene family, OsMATEl and OsMATE2 through expression in heterologous system, Arabidopsis. 
Expression of OsMATEs in Arabidopsis altered growth and morphology of transgenic plants. Genome-wide 
expression analysis revealed modulation of genes involved in plant growth, development and biotic stress in 
transgenic lines. Transgenic plants displayed sensitivity for biotic and abiotic stresses. Elevated pathogen 
susceptibility of transgenic lines was correlated with reduced expressions of defence related genes. Promoter 
and cellular localization studies suggest that both MATEs express in developing and reproductive organs 
and are plasma-membrane localised. Our results reveal that OsMATEl and OsMATE2 regulate plant 
growth and development as well as negatively affect disease resistance. 

Multidrug and Toxic compound Extrusion proteins (MATE) are a family of secondary active transporters 
which utilize electrochemical gradient of membrane maintained by ATPases for their transport activity 1 . 
Members of this transporter family have ubiquitous occurrence in all the organisms including human 
and plants. The first member of this family (NorM) was characterized from Vibriopara haemolyticus, which 
effluxes nofloxacin and ciprofloxacin outside of the cells in an energy dependent way 2 . Most of the earlier studies 
suggest that MATE proteins act as an efflux pump that export majority of drugs and xenobiotic compounds 
outside the cell, and largely contribute to the drug resistance in bacteria. Recent reports indicate that MATE 
proteins function in antiport manner and drive the substrate transport in exchange with H + /Na + depending on 
the organism. Studies suggest that H + coupling is operational in plants, which is substituted by Na + in bacteria. 
Though structure and mechanism of action of MATEs have not been studied in detail, recent investigations 
revealed conformation and transport behaviour of MATE proteins 3,4 . Conversely to designated name "multi- 
drug", numerous studies revealed that MATE proteins have stringent substrate specificity and facilitate the 
movement of specific compounds 5 . Presence of homologous putative MATEs has been identified from several 
organisms through comparative genomics. A genome wide scanning and analysis revealed that 58 MATE 
paralogues are present in Arabidopsis genome 6 . In general, MATE proteins from various source organisms have 
a conserved domain and share at least 40% amino acid sequence homology 3 . 

In higher plants, studies demonstrated that MATE proteins are mainly involved in the transport and trafficking 
of xenobiotic and small organic molecules 7 . The foremost characterized Arabidopsis MATE protein, AtDTXl, 
was demonstrated to export norfloxacin due to its ability to restore drug tolerance in the norfloxacin sensitive 
bacterial mutant (kam3 mutant) 8 . Recently, evidences for role of MATEs in the plants are rapidly accumulating. 
Depending on the established functions of this gene family in plants, MATEs have been grouped into three major 
classes. The first group of the MATE has been demonstrated to be associated with disease resistance in 
Arabidopsis. Several Arabidopsis mutants have been identified and examined for understanding the molecular 
mechanism of pathogen resistance. Among them, enhanced disease susceptibility mutant (eds5) of Arabidopsis 
displayed a reduced basal resistance during pathogen interaction. The fine genetic mapping of the eds5 locus 
identified that EDS5, a member of the MATE transporter family, was localised within this locus, and substantially 
contributed to disease tolerance 9 . The relation of EDS5 with disease tolerance has been explored very recently and 
demonstrated that EDS5 export salicylic acid (SA) outside chloroplast, where SA synthesis takes place 10,11 . 
Similarly, ADS1 (activated disease susceptabilityl), a negative regulator for the disease resistance, was another 
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MATE identified in Arabidopsis showing involvement during patho- 
gen infection 12 . A second group of MATE proteins is known for 
exporting small organic molecules such as citrate outside the cell 
which acts as a ligand molecule to bind aluminium (Al) in the rhizo- 
sphere. Unlike to other metals, Al toxicity is much pronounced in 
acidic soil and several genetic studies were performed to map the 
locus responsible for Al tolerance in plants. A locus for Al tolerance, 
Alt SB , was identified that encode a MATE via positional cloning in 
sorghum. This SbMATE has been shown to efflux the citrate from cell 
and forming non-toxic complexes with Al in soil solution 13 . Similar 
to sorghum, MATE proteins from Hordeum 14 , Arabidopsis 15 , 
Triticum 16 and maize 17 have been shown to participate in Al tol- 
erance. Tandem duplication of MATE paralogues in the maize gen- 
ome have been demonstrated to offer another level of transcriptional 
regulation to provide adaptation of maize to Al toxicity 18 . 

Another, third group of MATE transporters are involved in traf- 
ficking of secondary metabolites in vacuoles. Usually, vacuoles are 
the major repository site of most of the conjugated form of flavonoids 
comprising mainly flavonols, anthocyanins and flavone glycoside. 
The transport and storage of these flavonoides into vacuole are 
mediated by the different class of transporters on tonoplast 19 . A 
member of MATE gene family (TT12) was elucidated for sequest- 
ration of proanthocyanidins in vacuoles of seed that leads to pig- 
mentation of seed coat 20 . The kinetic study of vesicles isolated from 
TT12 expressing yeast revealed that it can specifically transport gly- 
cosidic form like epicatechin 3'-0-glucoside and cyanidin 3-0 gluc- 
oside 5 . Consistent to the TT12, MATEs from Medicago and Vitis 
have been characterized for mobilization of flavonoids in cell orga- 
nelles 21 23 . Similar to the role in vacuolar flavonoids sequestration, 
MATE proteins are also recognised for alkaloids trafficking into 
vacuoles in tobacco 24,25 . 

Apart from these functions FRD3, a MATE citrate efflux trans- 
porter, was shown to involved in Fe transport towards stele region 
and hence assisted in distribution of Fe throughout aerial parts in 
Arabidopsis 26 . Curiously, one member of MATE gene family (Mdtl) 
was shown to regulate the glucose assimilation, and which is thought 
to be necessary for nutrient utilization, sporulation and pathogeni- 
city in fungi 27 . In an another report, role of a MATE-like transporter 
(ZRZ) has been demonstrated in organ initiation as overexpression 
of this MATE led to diverse morphological changes in Arabidopsis 26 . 
Therefore, it seems that MATEs govern the diverse physiological 
functions in the plants which need to be elucidated. Likewise to 
Arabidopsis, genome-wide analysis of MATEs encoding loci in rice 
revealed the presence of large members of this gene family however, 
function of only a few members has been elucidated as yet 2 " 30 . Our 
earlier study suggested that expression of some of the genes encoding 
MATEs is up-regulated in rice seedlings exposed to arsenic (As) 
stress 31 . Of these loci, we have investigated the role of two gene, 
OsMATEl (Os03g08900) and OsMATE2 (Os05g48040) using het- 
erologous system, Arabidopsis, in this study. Our observations sug- 
gest that these MATEs play crucial role in plant growth, development 
and in stress responses. 

Results 

In silico identification and analysis of rice MATEs. In order to 
identify MATE genes present in the rice genome, BLASTP analysis 
at Rice Genome Annotation Project using Arabidopsis MATEs as 
query sequences was performed. A total of 53 members of the rice 
MATE family were identified (Supplementary Table S2). To study 
the evolutionary relationship between identified rice MATE pro- 
teins, a maximum likelihood phylogenetic tree was constructed 
using all the MATE proteins from Arabidopsis as well as from 
other plants with well established functions. The analysis suggests 
that rice MATE proteins clustered in different groups of Arabidopsis 
and other known MATEs from different plants. 



As very limited information about the function of rice MATE 
proteins is available, we have selected two arsenic responsive mem- 
bers 31 , OsMATEl (Os03g08900) and OsMATE2 (Os05g48040), for 
functional studies. Though both of MATEs show same responses 
with As, whereas grouped in two distinct clades (Supplementary 
Fig. SI). These rice MATEs, OsMATEl and OsMATE2 encode poly- 
peptide comprising the 489 and 500 amino acid residues respectively 
and shares 36% identities and 53% similarity with each other within 
90% query cover. Furthermore, trans-membrane helix prediction 
through TMHMM indicated the presence of 12 trans-membrane 
helices in both MATEs (Supplementary Fig. S2). This preliminary 
analysis suggests that OsMATEl and OsMATE2 were membrane 
bound proteins and possibly engage in transport activity across 
membrane. 

Expression of OsMATEl and OsMATE2 alter growth in Arabi- 
dopsis. In order to investigate the function of OsMATEl and 
OsMATE2, full-length cDNAs were expressed in Arabidopsis under 
control of a CaMV35S promoter. The presence and expression of 
transgene was confirmed through genomic and semiquantitative RT- 
PCR (Supplementary Fig. S3). Depending upon relative expression 
and phenotypic appearance, three independent lines expressing 
OsMATEl and OsMATE2 were selected for further study. Up to 
two weeks of germination in soil, growth pattern of the all trans- 
genic lines was similar to WT plants. However, after two weeks of 
growth, transgenic lines expressing OsMATEl or OsMATE2 displayed 
a distinguishable growth pattern in comparison to WT plants, but 
very similar to each other. The visible changes at the stage of three 
week old plants among all transgenic lines were quite prevalent 
(Supplementary Fig. S3, Fig. 1) that consist petiole length, leaf size 
and rosette leaf arrangement in vegetative parts. 

The pattern of rosette leaves of the mature plants was different in 
comparison to WT plants without any change in the number of 
leaves (Fig. la). Petiole length and rosette leaf size of the transgenic 
lines were significantly increased compared to the WT (Fig. lb, c). 
Collectively, the expression of OsMATEl and OsMATE2 resulted 
variation in vegetative tissues, including the leaf size and altered 
pattern of rosette development. To authenticate these inferences, 
growth of two T-DNA insertional mutants of Arabidopsis MATEs 
(SALK_045655C and SALK_124549C) having close sequence sim- 
ilarity with OsMATEl and OsMATE2 were analysed. Interestingly, a 
cognate growth pattern of salk lines with WT plants was observed 
(Supplementary Fig. S4a). Altogether, these observations suggest that 
OsMATEl and OsMATE2 play important role in the growth and 
development of plants. 

Altered floral morphology in transgenic lines expressing OsMATEl 
and OsMATE2. Aside to morphological changes, variations in flower- 
ing time and flower size were observed in all the transgenic lines 
expressing OsMATEl and OsMATE2 in comparison to WT plants 
(Fig. 2). An early flowering was observed in transgenic lines under 
normal growth conditions. Flowering time in transgenic lines was 
significantly reduced (up to four days) relative to WT plants (Fig. 2) 
but rosette leaf numbers at bolting stage were identical in transgenic 
lines and WT plants. Flowering time of salk mutants and WT plants 
were found similar and no significant difference was observed during 
initiation of bolting (Supplementary Fig. S4b). In addition to the 
alteration in flowering time, significantly reduced size of mature 
flower (length and width) was observed in transgenic lines com- 
pared to WT (Fig. 3). Though silique size of each transgenic line 
was comparable to WT, the number of siliques per plant were 
significantly higher (up to 50%) in comparison to WT (Fig. 4b). 
Transgenic lines also produced smaller seeds as compared to WT, 
although the viability of seeds was not compromised (Fig. 4b). 
Similar to these differences of reproductive organs, seed number 
produced per silique in all the transgenic lines was significantly 
higher corresponding to WT plants (Fig. 4c). Taken together, these 
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Figure 1 | Alteration of leaf morphology in transgenic lines. (A) Leaf sizes and pattern of each rosette leaf beginning from the oldest one of WT and one 
line representing OsMATEl and OsMATE2 transgenic plants. Plants were grown for 4 week under 16 h light-8 h dark photoperiod before photograph. 
(B) Photograph showing the difference in leaf size (9" 1 leaf of rosette from oldest one) of WT and 3 representative lines of OsMATEl and OsMATE2 
transgenic plant. White line indicates the scale bar (1 cm). (C) Ratio of length and width of leaf (9 th position) of transgenic lines and WT. Bars are mean 
and error bars represent ± SE of 8-10 individual plants. Level of significance was evaluated by One-way ANOVA (Newman-Keuls used as post hoc test). 
***Significantly different from WT (P < 0.0001). 



results suggest that expression of OsMATEl or OsMATE2 in Arabi- 
dopsis develops many phenotypic changes which include reduced time 
of transition from juvenile to flowering stage, morphology of flowers, 
silique number, seed size and silique number produced per plants. 
However, such a macroscopic changes were not observed in salk lines. 

Genome wide expression analysis in the transgenic lines. To add- 
ress reason for phenotypic changes in transgenic plants, modulation 
of transcriptome was analyzed using Affymetrix ATH1 (22K) array 
of three week old WT and transgenic line expressing OsMATEl. At 
this stage, clear phenotypic changes between WT and transgenic 
lines were visible. After normalization of eel files through d-Chip, 
a significantly modulated expression (1.6 2 fold, p values 0.05) of 39 
genes in transgenic lines compared to WT plants were observed. Of 
these, expression of 14 and 25 genes were up- and down-regulated, 
respectively in the transgenic line (Table 1 and 2). The majority of 
genes up-regulated were related to DNA binding proteins and 
transcription factors controlling the diverse metabolic and 
physiological processes in plants. Out of most notable modulated 
expression, pentatricopeptide (PPR) repeat-containing protein up- 
regulated nearly 3 fold (Table 1). Characterized PPR proteins in 
Arabidopsis are known to be localized in mitochondria and 
participate in the post transcriptional modification of RNA. The 
impaired or loss of function of PPR lead to reduced NADH 
dehydrogenase activity and consequently a compromise defence 
response in Arabidopsis 32,3 *. Enhanced expression of RD20 (Respon- 
sive to Desiccation 20) and CAX7 (calcium exchanger 7) was also 



seen, which are involved in salinity stress. A significant up-regulation 
of NAC6 and NAC32 was observed in transgenic plants. NAC 
domain containing proteins are highly studied transcription factor 
family in plants and are known to regulate growth and development 
of plants and drought, desiccation, salinity and ABA response 34 . 
Interestingly, expression of a Jumonji transcription factor (jmjC 
domain-containing protein) was enhanced in transgenic lines. In 
Arabidopsis, jmjC domain-containing protein (Jmj30) is known as 
the key determinant of circadian rhythm and regulates the flowering 
time in a cooperative manner with central oscillator, circadian clock 
associatedl and late elongated hypocotyl 35 . Expression of all of the 
up-regulated genes was further validated through qRT-PCR (Fig. 5) 
that followed the same pattern of expression as observed in micro- 
array analysis. Overall, these observations revealed modulated expres- 
sion of certain regulators involved in contending the abiotic stress, 
and one specific regulator of circadian function (Jmj30) in transgenic 
plants. 

Our analysis also indicated that a group of defence related genes 
were down-regulated in transgenic plants in comparison to WT 
(Table 2). Quantitative RT-PCR analyses of these genes were in 
agreement with microarray results (Fig. 5). These genes comprised 
the pathogenesis related proteins (PR-1, PR-5), R-l,3-glucanase, 
chitinase, plant natriuretic peptide A, aspartyl protease family pro- 
tein, crinkly related protein and ankyrin binding protein. A 
remarked decrease in transcript level of three cysteine rich receptor 
like protein kinases (CRK) and two receptors like protein were 
observed among 28 down-regulated genes. PR proteins are, in 
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Figure 2 | Early flowering in transgenics lines expressing OsMATEs. (A) Expression of either of OsMATEl or OsMATE2 in Arabidopsis induces early 
flowering in Arabidopsis. Photograph of two lines of OsMATEl and OsMATE2 (nearly 3 week old) along with WT are displayed. (B) Number of rosette 
leaf and days at bolting of WT and 3 transgenic lines of the OsMATEl and OsMA TE2 grown under LD conditions (16 hlight-8 h dark cycle). Data are the 
mean and ±SE (n = 8-12). Significance was evaluated by One-way ANOVA (Newman-Keuls used as post hoc test). * "Significantly different from 
WT (P < 0.0001). 



general, recruited downstream in salicylic acid signalling pathway 
and extend the cascading of signals 9 . Expression of (3- 1 ,3-ghicanase, a 
class of PR proteins, and chitinase are induced during fungal infec- 
tion that integrate a sort of tolerance while invasion of fungus in the 
plants 12 . We conceived that the significantly reduced expression of 
stress related genes may lead to compromised defence response in 
transgenic plants. 

Arsenic response of transgenic lines. Our earlier study suggested 
that OsMATEl and OsMATE2 expression is enhanced during As 
exposure in rice. To examine the response of As on growth of 
OsMATEs expressing transgenic lines, seeds of transgenic lines, 
WT and pBI121 (empty vector transformed control) were grown 
for ten days on half strength MS containing agar plates with and 
without addition of As. All the genotypes grew similarly in absence 
of As(V). However, growth and root length of transgenic lines were 
significantly decreased in comparison to WT and pBI121 plants 
under As(V) rich environment (Supplementary Fig.S5). As partici- 
pation of MATE proteins in combating the aluminium ( Al) stress is 
well documented, the effects of Al on growth of transgenic lines were 
also evaluated. However, root length of all the genotypes was same at 
presence of Al in media and no change has been observed compared 
to WT and pBI121 plants (Supplementary Fig. S6). The root length 
data suggests that As(V) specifically sensitizes the growth of 
transgenic lines. 

Enhanced disease susceptibility in transgenic lines. Expression 
analyses suggested that a group of defence related genes were 



down-regulated in transgenic plants compared to WT (Table 2). 
We envisioned that significantly reduced expression of stress 
related genes may lead to compromised defence response in trans- 
genic plants. To further study the expression of MATEs under biotic 
stresses, we have investigated the expression pattern of OsMATEl 
and OsMATE2 in Nipponbare and IAC165 cultivar of rice during 
infection of M. grisea isolate (FR13) and S. hermonthica using 
publically available microarray data 36,37 . Expression dataset exhibited 
that OsMATEl was enhanced at post infection of M. grisea whereas its 
expression was repressed or unaffected during S. hermonthica interac- 
tion. By contrast, OsMATE2 induced specifically by the exposure of 
S. hermonthica (Fig. 6a). The expression analysis revealed that 
OsMATEl and OsMATEl were differentially expressed during M. 
grisea and S. hermonthica infections in rice and suggest that these 
genes might have a role in biotic stresses. 

To ascertain the role of these MATEs in biotic stress, basal resist- 
ance of the transgenic and WT plants challenged with Pseudomonas 
syringeae (PstDC3000), a natural host of Arabidopsis causing spek 
disease was measured. After 3 dpi, an explicit enhanced susceptibil- 
ity was observed in transgenic lines in comparison to WT in the 
form of the chlorotic lesion with a substantially higher bacterial titre 
(Fig. 6b, c). Our observations suggest that OsMATEl and OsMATE2 
may be negatively regulating disease tolerance apart from governing 
the growth and development. Furthermore, we have carried out 
infection assay after sprayed inoculation of bacterial culture on whole 
plants and quantified the bacterial titre at 7 dpi, which reconfirmed 
the susceptibility and significantly higher titre in all the transgenic 
lines expressing OsMATEl and OsMATE2 (Supplementary Fig. S7). 
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Figure 3 | Variation of flower size in transgenic lines. (A) Flower size, 
both length and width, was reduced in transgenic plants. Photograph of 
flowers were taken after 5 week growth under standard conditions. (B) Bars 
are mean of length and width of 3 independent lines of each gene and 
variation represent as ±SE (n = 8-12). Significance was evaluated by One- 
way ANOVA (Newman-Keuls used as post hoc test). ***Significantly 
different from WT (P < 0.0001). 

Promoter of OsMATEs is active in seedling and flowering tissue. 

Expression of OsMATEl and OsMATE2 in Arabidopsis led to several 
alterations related to plant growth and development. To establish 
function of these genes in detail, we studied tissue-specific expres- 
sion of uidA gene under control of OsMATEl and OsMATE2 
promoters in Arabidopsis. In this regard, in silico analysis of 
promoters (ProMATEl, 1540 bp; ProMATE2, 1434 bp) has been 
carried out, which indicated the presence of various conserved 
consensus sequences of ris-regulatory elements. ProMATEl has the 
41 unique motifs which were absent in ProMATE2, Similarly, 
ProMATE2 contains 30 unique motifs (Fig. 7a). However, 36 
common motifs were found in promoter region of both the genes 
(Fig. 7a). The detailed descriptions of these motifs are provided in 
Table S3, S4 and S5. Examination of GUS activity in all the tissues 
representing various developmental stages of Arabidopsis revealed a 
prominent GUS staining in newly formed leaves, shoot apex and in 
the root of ten days old seedlings for both OsMATEl and OsMATE2 
promoters (Fig. 7b). More importantly, no GUS activity was 
observed in other studied tissues (rosette and cauline leaves, stem) 
except in the reproductive tissues. After the fertilization, a strong 
GUS activity was seen in immature flower buds, early stage of 
embryo and seed development (Fig. 7b). Taken together, histolo- 
gical GUS staining of OsMATEl and OsMATE2 promoters in 
Arabidopsis showed a cognate expression pattern, and suggests 
that reporter activity may be driven by the common motifs present 
in these promoters. Furthermore, a strong promoter activity at the 
shoot apex and during reproductive phase suggests the likely role of 
OsMATEl and OsMATE2 in plant growth and development. 

Sub-cellular localization of OsMATEl and OsMATE2. In order to 
investigate the sub-cellular localization of OsMATEl and OsMATE2, 
complete coding region of OsMATEl and OsMATE2 were cloned in 



transient expression vector 326 sGFP 38 and transformed in protoplast 
isolated from mesophyll cells of Arabidopsis leaves. The examination 
of GFP fluorescence, after overnight incubation in dark, under 
confocal microscope explicitly suggests that OsMATEl and OsMA- 
TE2 underlie on the plasma membrane, whereas a diffused green 
fluorescence was seen throughout the cell transformed with empty 
vector (Fig. 8). Altogether, these observations indicated that 
OsMATEl and OsMATE2 reside on the plasma membrane. 

Altered flavonoid content in transgenic lines. The relatedness of 
OsMATEl and OsMATE2 with some of characterized flavonoid 
transporters in phylogenetic tree (Supplementary Fig. SI) suggest 
the possible role of these proteins in flavonoid transport. The two 
fold higher expression of GST (Glutathione-S transferase) in 
OsMATEl expressing lines as observed in expression analyses also 
strengthen this hypothesis. The physical conjugation of glutathione 
to flavonoid catalysed by GST is a crucial metabolic step for 
flavonoids trafficking 39 41 . Therefore, content of major flavonols 
(kaempferol, quericetin and rutin) have been quantified in leaves 
of WT and transgenic lines. Flavonol profiling revealed the perturba- 
tion of kaempferol content in transgenic lines in comparison to WT 
plants. As speculated, kaempferol content was enhanced in 
transgenic leaves, which was undetectable in WT leaves (Table S6), 
and suggests that a proper flavonol channelling might be affected in 
transgenic plants. 

Discussion 

By using reverse genetics approach and a combination of molecular 
and biochemical tools, functional characterization of two members 
of the MATE family from rice, OsMATEl and OsMATE2 has been 
carried out in this study. The scarcity of knowledge about the func- 
tion of MATE proteins in rice, except the known citrate efflux activity 
of the FRDs 29,30 , encouraged us to ascribe the detailed role of 
OsMATEs out of total 53 proteins in rice genome. The constitutive 
expression of two As-responsive MATE encoding genes, OsMATEl 
and OsMATE2, led to the diverse changes in growth and morphology 
of Arabidopsis. The major alterations were seen in leaf morphology, 
pattern of rosette leaves arrangement and flower development. We 
further demonstrate that transgenic lines were sensitive to As(V) 
exposure as well as exhibit a compromised basal resistance. 
Repressed immunity of transgenic lines was correlated with lowered 
expression of genes encoding pathogenesis related proteins. 

Astonishingly, transgenic plants harbouring OsMATEl and 
OsMATE2 independently displayed distinct common alterations 
that encompass leaf size, rosette leaf arrangement, flower, silique 
and seed morphology (Fig. 1-4). Nonetheless, these phenotypic var- 
iations were lacking in salk lines possibly due to presence of large 
members of this gene family with overlapping functions in 
Arabidopsis. For studying the mechanisms and responsible genes, 
we studied the genome wide expression analyses of transgenic and 
WT leaves. The analysis suggested that enhanced expression of a set 
of genes in transgenic plants may play important role in such pheno- 
typic changes (Table 1 and 2). Two loci encoding the protein that 
contains NAC domain representing NAC transcription factor family 
were unregulated. NAC proteins constitute one of the largest families 
of transcription factors in plants and NAC domain containing tran- 
scription factors is considered as key regulator of stress perception 
and co-ordinating developmental programmes in plants 42 . Recently 
phytochrome-b interacting proteins, vascular plant one zinc fingerl 
(VOZ1) and VOZ2, were identified to regulate the flowering in 
Arabidopsis. VOZ1 and VOZ2, a transcriptional activator, belong 
to the NAC family and vozlvoz2 double mutant showed a late flower- 
ing phenotype 43 . In this study, we also reported an early flowering in 
transgenic lines and changes in silique number, flower and seed size 
corresponding to WT plants (Fig. 2-4). These changes might be due 
to increased transcript level of NAC transcription factors. Early 
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Figure 4 | Alteration in silique and seed number in transgenic plants. (A) All the transgenic lines more silique per plant than WT. Photograph of 6 week 
grown plants (2 lines of each transgenic) was taken. Data of 3 transgenic lines presented as mean and ±SE (n = 8-12). (B) Seed size of transgenic and WT 
plants (representative of 1 line of each transgenic). (C) Seeds number per silique was high in transgenic lines. Values were collected from 10 silique 
per plant and 8-12 plant for each lines and WT. Data are presented as mean, ±SE (n = 8-12). Significance was determined by One-way ANOVA 
(Newman-Keuls used as post hoc test). ***Significantly different from WT (P < 0.0001). 



flowering of transgenic plants can be due to enhanced expression of 
JumonjiC (JmjC) containing protein. Basically, Jmj30 is a type of 
histone demethylase, which epigenetically regulated the circadian 
clock of Arabidopsis by affecting histone remodelling 35 . In 



Arabidopsis, circadian clock regulates various processes such as gene 
expression, photoperiodic flowering, and leaf movement 44 . Thus, 
higher expression of Jmj30 might have a role in association with 
other factors as a determinant of early flowering in transgenic plants. 



Table 1 | List of 
Locus ID 


significantly up-regulated genes in OsMATEl expressing Arabidopsis lines 

Annotation 


Fold change 


At5g593 1 0 


LIPID TRANSFER PROTEIN 4 


6.87 


At5g39610 


NAC domain containing protein 6 


3.02 


At3g20810 


JUMONJI (jmjC) domain-containing protein 


2.79 


Atlg07590 


PENTATRICOPEPTIDE (ppr) repeat-containing protein 


2.75 


At2g 19970 


Pathogenesis related protein, putative 


2.52 


At5g 17860 


CAX7 (calcium exchanger 7) 


2.24 


At2g33380 


RD20 (response to dessication20) 


2.19 


Atlg30260 


Unknown protein 


2.18 


AT2G29490 


Glutathione S-transferase TAU 1 


2.03 


Atlg77450 


NAC domain containing protein 32 


2.02 


At5g42900 


Unknown protein 


2.02 


At3g28290 


UPF0496 protein 


2 


At2g36120 


Glycine-rich_protein 


1.99 


At4g 15490 


UDP-glycosyltransferase 


1.97 


Atlg77760 


Nitrate Reductase 1 


1.95 


At5g61820 


Unknown protein 


1.95 


Atlg72260 


Toxin receptor binding (thionin2.1) 


1.83 


At5g23240 


DNAJ heat shock N-terminal domain-containing protein 


1.68 
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Table 2 List of signi 
Locus ID 


■icantly down-regulated genes in OsMATEl expressing Arabidopsis line 

Annotation 
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Plont Ccl rssistoncs protsin 
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At9n 1 RAAO 
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pLnf nntri 1 1 rati/** nantirJa A /Ph\IP Al 
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—4 ?4 


At^nl 07A0 
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-4 2 1 
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Po/-ar\t.-ir 1 iba Pr/-\toin 9 T /PI P9T\ 

Kecspror liks rrorein zo [KLrzoj 
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Atl n?0900 


VACIJOIAR SORTING RFrFPTOR A (VSRA1 


-2 77 
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rtlNIx lUlllvyllill piUlclll UlilUlliy 


-9 69 
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Atlg35230 
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Cysteine-rich receptor-like protein kinase 4 


-2.02 


At5g24530 


Flavanone 3-hydroxylase-like protein 


-1.99 


Atlg35210 


Unknown protein 


-1.97 


At3g27690 


Photosystem II light harvesting complex gene (LHCB2.3) 


-1.9 


At2g32030 


GCN5-related N-acetyltransferase (GNAT) family protein 


-1.81 



We further checked the co-expressed genes of Jmj30 by searching at 
ATTED-II to get a more clue at molecular level. Our analysis suggests 
that C0L9 ( At3g07650) and ELF4 ( At2g40080) show a high correla- 
tion co-efficient along with Jmj30. C0L9 belongs to CONSTANS 
gene family and is a negative regulator of photoperiodism and 
flowering via suppressing the expression of Flowering locus T (FT) 
and Suppressor of Overexpression of C01 (SOC1) 45 . A subtle 
increase in expression of C0L9 in transgenic lines supports the 
assumption (Supplementary Fig. S8). In addition, induced express- 
ion of the pentatricopeptide repeat containing (PPR) and other redox 
related genes suggest an imbalance of proper redox state of the cells. 



There are few recent report emerges that correlate a synchrony 
between oxidative burst (ROS signalling) and circadian clock of 
Arabidopsis 46 . 

Our expression analyses revealed that decreased expression of a 
group of defence related genes in OsMATEl expressing Arabidopsis 
lines compared to WT (Table 2). These genes can be grouped in two 
broad categories, some of genes contribute to restrict the invasion of 
pathogen as their expression was highly induced, mainly comprising 
the PR proteins, (3-1,3-Glucanase and Receptor like proteins 
(Table 2). Proteins encoded by these genes play critical role in patho- 
gen recognition and signal perception followed by activation of 
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Figure 5 | Modulation of gene expression in transgenic line expressing OsMATEl. Relative gene expression between WT and transgenic plant was 
calculated using qRT-PCR ( tubulin was taken as endogenous control). Data provided represent fold expression Log 2 ( transgenic/ WT ) . Specific-primer set 
of individual genes was designed from intron spanning region. 
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Figure 6 | Bacterial infection assay. (A) Expression of OsMATEl and OsMATE2 during infection of M. grisea on Nipponbare and S. hermonthica on 
Nipponbare and IAC165 cultivars of rice (Oryza sativa L.) at various indicated time, M2-M11 (mock) and D2-D11 (infection), days post inoculation of 
Affymetrix Rice Genome Array as cluster display. The colour scale (representing log signal values) is shown at the bottom. (B) Symptom of disease in 
leaves of WT and two transgenic lines at 3 dpi after PsfDC3000 spot inoculation. (C) Bacterial titre in the leaves was quantified at 3 dpi as shown as mean 
and standard errors calculated from three biological replicate (each leaf from different plants). The experiment was carried out more than three times. 
Significance was determined by One-way AN OVA (Newman- Keuls used as post hoc test). *** indicates significant difference (P < 0.0001) as compared 
to WT. 



defence responses. Similarly, cysteine-rich receptor like kinases and 
leucine rich repeat (LRR) protein kinases were significantly down- 
regulated (Fig. 5). A NB-LRR is the major class of plant receptor, 
which intercepts effectors signal from much of pathogens and con- 
trol the transcriptional mobilization of defence pathways in 
Arabidopsis 47 . Reduction in the expression of various genes respons- 
ible for signal perception to many downstream processes reflect the 
reduced basal disease resistance in transgenic plants and that was 
manifested in PsfDC3000 infection assay as well (Fig. 6). An 
enhanced susceptibility was observed in transgenic lines expressing 
OsMATEl and OsMATE2 compared to WT plants either on the 
detached leaf (Fig. 6b, c) or after seven dpi in planta sprays of bac- 
terial inoculums (Supplementary Fig. S7). Similar to our observa- 
tions, another study revealed that an activation tagging mutant of 
adsl-D (MATE protein) is highly susceptible to PsfDC3000 and 
knockdown of ADS1 results a higher basal resistance in 
Arabidopsis 12 . A systematic proteomic analysis of root exudates of 
Arabidopsis revealed that chitinase, (3-1,3-Glucanase and other 
defence proteins are secreted in higher amount during flowering 
and shown to link between flowering and defence processes 48 . The 
involvement of MATE proteins in disease resistance is unexception- 
ally accepted, and formerly identified during genetic screen of 
Arabidopsis mutants 912 however underpinning molecular mech- 
anism is explored very recently 10-11 . It has been revealed that EDS5 
regulate salicylic acid biosynthesis by exporting from chloroplast. 
Flavonoids are often regarded as prominent signalling compound 
and interplay during biotic stresses however underlying actual 
molecular events of flavonoids signalling is not very much known. 
It has been reported that through managing accumulation of some 
flavonols plants can better adapt under plenty of pathogens and 
insect attacks 49,50 . In our study, flavonoid profiling of transgenic 



and WT leaves suggested altered flavonol content in transgenic 
plants in comparison to WT plants. We propose that perturbation 
of kaempferol, through OsMATEs expression in Arabidopsis, could 
be the central and foremost event that persuades morphological 
alterations. Although, specific substrates transported by these trans- 
porters have not resolved yet, but the results presented here indicated 
that it may be a flavonol, and which need to explore in future by 
molecular transport experimentation. 

A prominent (AsV) sensitivity was observed in transgenic lines 
expressing MATEs as compared to WT plants (Supplementary Fig. 
S5). The sensitivity of transgenic plants could be explained on the 
basis of inevitable contribution of flavonoids in maintaining redox by 
retaining ROS scavenging activity. Disturbance of this process via 
perturbation of intracellular flavonoids either at different steps of 
their biosynthesis or transport could lead to the manifestation of 
As(V) sensitive phenotype. Plasma- membrane localized MATE 
transporters are devoted for citrate efflux while others are recognised 
for flavonoid and other secondary metabolites sequestration in 
vacuoles 25,30,51 . The previously identified role of MATE transporters 
engage in flavonoids sequesteration are found to mostly restricted on 
vacuolar membrane rather than plasma-membrane localization of 
OsMATEl and OsMATE2, which is elucidated in this study (Fig. 8). 
Thus, more and more information generated about functional sig- 
nificance of members of this gene family in plants, would elaborate 
our insight about transport of flavonoids and other compounds ran- 
ging from intracellular to long distance transport and their impact on 
various physiological processes in future. 

In conclusion, our study suggests that OsMATEl and OsMATE2 
are functionally redundant to each other based on observations 
obtained in transgenic lines as both MATEs direct the common 
phenotypic changes through expression in Arabidopsis as well as 
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Figure 7 | Promoter analyses of OsMATEl and OsMATE2. (A) cis-regulatory elements of ProOsMATEl and ProOsMATE2 was analysed through 
PLACE and venn diagram shows the number of common and unique motifs present in promoters. (B) Promoter activity in WT and Arabidopsis lines 
expressing ProOsMATEl :uidA and ProOsMATE2:uidA. Representatives of GUS staining performed in seedling (10 day old), flowers and silique of mature 
plant. White scale bars present in each panel in (B) are equivalent to 100 urn length. 



their similar promoter activities in Arabidopsis. There is a necessity 
of thorough study in the future to obtain more insight about nature of 
substrates transported by these proteins to establish a direct link 
between growth-development and disease resistance in plants. 

Methods 

Plant material and growth condition. Arabidopsis thaliana (L.) accession Col-0 and 
salk mutants (SALK_045655C and SALK_124549C) were used in the study. 
Arabidopsis plants were grown in growth chamber (Canviron, USA) under controlled 
conditions 16 h light-8 h dark photoperiod cycle, 22°C temperature, 150-180 umE 
m~ 2 s -1 light intensity and 60% relative humidity. Seeds were sown in soilrite and 
stratified for 48 h in the dark at 4 C before shifting in the light. Plants were watered 
with nutrient media (5 mM KN0 3 , 2 mM MgS0 4 , 2 mM Ca (N0 3 ) 2 , 2 mM K 
H 2 P0 4 , 50 uM Fe(II). EDTA, 70 uM H 3 B0 3 , 14 uM MnCl 2 , 1 uM ZnS0 4 , 0.5 uM 
CuS0 4> 0.2 uM Na 2 Mo0 4 , 10 uM NaCl, 0.01 uM CoCl 2 ) and reverse osmosis (RO) 
water at alternate interval in every week until plants complete their lifecycle. 

Database search and Bioinformatics analysis. MATE sequences were searched from 
the database on the BLAST server (http://www.ncbi.nlm.nih.gov/blast/Blast.cgi and 
http://rice.plantbiology.msu.edu) by using keyword, domain name and BLAST 
against characterized Arabidopsis MATEs. We retrieved full length MATE sequences 
of rice, Arabidopsis, and fully characterized MATE proteins from other organisms. 
Multiple amino acid sequences were aligned with the help of ClustalW alignment tool 
(http://www.ebi.ac.uk/Tools/clustalw2/index.html). The phylogenetic tree using full- 
length MATEs was constructed by the neighbour -joining method using Phylip 52 and 
displayed using PhyloDraw V0.82 53 . Expression study of OsMATEl and OsMATE2 
during biotic stress in rice has been carried out using publicly available microarray 
CEL files in the GEO database. GEO database for virulent infection by Magnaporthe 
grisea (GSE7256) and expression data for interaction with the plant parasite Striga 
hermonthica (GSE10373) were used in expression analysis. The Affymetrix CEL files 
of each experiment were imported and analysed using d-Chip under default setting 54 . 

Cloning, construct preparation and transgenic plant development. Full-length 
cDNAs of OsMATEl and OsMATE2 were amplified from the single stranded cDNA 



library prepared from rice root though PCR by using nested primers (Table Si) and 
cloned in pT/Z subcloning vector (Fermentas, USA). Another set of primers was 
designed with restriction sites ofXbal and Sad enzymes in forward and reverse 
primers, respectively. The coding region of both MATEs was cloned in plant 
expression vector pBI121 between Xbal and Sad restriction sites. These constructs 
were transformed in Agrobacterium tumefaciens strain GV3101 by the freeze thaw 
method. Arabidopsis transgenic plants were generated through floral dip method 55 . 
Transformants were screened on half MS agar plates containing kanamycin (50 ug 
per ml media). Based on the relative expression of transgene through semi 
quantitative RT-PCR, three lines for each of OsMATEl and OsMATE2 among the 15 
screened independent lines were selected for further analysis. All the studies were 
performed in T4 homozygous lines. For promoter analysis, upstream sequences from 
the translation initiation codon of OsMATEl (1540 bp) and OsMATE2 (1434 bp) 
were retrieved and amplified using High Fidelity enzyme mix (Fermentas, USA) from 
genomic DNA and cloned in pCambial303 plant binary vector between Hindlll and 
Ncol restriction sites. Transgenic plants were generated and screened on hygromycin 
(10 jig per ml). Similarly, for preparing the construct employed in sub-cellular 
localization, coding region of OsMATEl and OsMATE2 were cloned in 326 sGFP 
transient expression vector between Xbal and BamHl restriction sites to encode the 
fusion proteins 56 . All the clones were sequenced to ensure the in frame fusion and 
error free preparation of constructs. 

Microarray and gene expression analysis. Microarray analysis was performed using 
RNA isolated from rosette leaves of nearly three weeks old Col-0 and OsMATEl 
expressing plants. Whole genome Arabidopsis ATH1 array (Affymetrix, USA) was 
utilized in experiment and each step was followed as per manufacturer's instruction. 
Total RNA was isolated from rosette leaves by using Spectrum Plant Total RNA Kit 
(Sigma, USA), followed by treatment with RNAase-free DNase (Fermentas, USA) to 
remove the genomic DNA. Briefly, 250 ng total RNA of each sample was taken as 
starting material and strictly followed the Affymetrix protocol during RNA 
amplification, hybridization, washing- staining and scanning of the arrays. 
Normalization of cell files and relative fold change were measured using d-Chip 54 . 
Quantitative Real time PCR (qRT-PCR) was carried out to validate microarray data. 
Total RNA (2 ug) was used in first strand cDNA library preparation (RevertAid H 
minus first strand cDNA synthesis Kit Fermentas, USA). One ul template was taken 
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Figure 8 | OsMATEl and OsMATE2 reside on plasma- membrane. Arabidopsis mesophyll protoplast transformed with 326sGFP (A-D), 
OsMATEl:GFP (E-H) and OsMATE2:GFP (I-L). Confocal cross sections show GFP fluorescence (excitation 488 nm and emission light measured 
through 505-530 nm band pass filter), autofluorescence of cell (excitation 543 nm, emission measured between 580 nm and 705 nm), bright field images 
and their overlay in different panels. (D) (H) and (I) represent overlay of (A-C), (E-G) and (I-K), respectively. Scale bars are included in images. 



in each reaction after 10 times dilution with nuclease free water. A reaction was setup 
with template, forward primer, reverse primer and Fast syber green mix in AB Fast 
7500 thermal cycler (Applied Biosystems, UK). Unique primer set (MWG, India) was 
designed from the intron spanning region of each studied gene, wherever possible, 
and tubulin was taken as endogenous control to normalize the data (Supplementary 
Table Si). For semi -quantitative RT PCR analysis, first strand cDNA was taken as a 
template and a primer set specific for OsMATEl and OsMATE2 (Table SI) and 2X 
master mix (Fermentas, USA) were used. A reaction set up, with an initial 
denaturation at 94 C for 2 min, 28 cycles at 94 C for 15 s, 55 C C for 15 s, and 72' C for 
30 s, followed by a final 5 min extension at 72°C, condition were executed in 20 ul 
reaction volume. PCR products were separated in 1.2% agarose gel in 0.5 X TBE 
running buffer. For taking equal amount of cDNA tubulin was taken as an 
endogenous control. 

Metal exposure and root growth assay. Seeds of all the genotypes were surface 
sterilized with sodium hypochlorite and plated on half strength MS agar plates with 
and without supplementation of As(V) and Al 57 . Plates were kept in dark for 48 h at 
4°C (stratification) and photographed after the growth of ten days in standard growth 
conditions. Root lengths of each genotype and WT plants were manually measured 
and taken from at least three separate plates. 

Analyses of transgenic lines. Homogygous transgenic lines expressing OsMATEl 
and OsMATE2 were grown for complete life cycle and all observations regarding 
different phenotypic characters were collected at indicated growth stages and 
monitored at least 2-3 successive generations. For histological GUS activity, 
ProOsMATET.uidA and ProOsMATE2:uidA transgenic lines were grown on half 
strength MS agar plates for 10 days or on soilrite. Seedlings and different tissues from 
mature plants were rinsed in ice cold PBS buffer 3-4 times and after three washing 
samples were incubated 16 h at 37°C in GUS staining buffer containing 100 mM 
NaP0 4 , 1.0 mM ferrocyanide, 1.0 mM ferricyanide, 0.05% TritonX-100,pH7.2, and 
2 mM 5-bromo-4-chloro-3-indolyl-b-D-glucuronide (initially dissolved in N,N- 
Dimethyl formamide) 58 . After successive washing from the gradient of ethanol (70%- 
90%), images were captured under light microscope (Carl Zeis, USA). 

Bacterial infection assay. Pseudomonas syringae pv. tomato strain (PsfDC3000) was 
allowed to grow overnight in king's B media containing 1 5 ug rifampicin per ml 
media at 28"C with 220 rpm shaking. Cells were collected by 5000 rpm 
centrifugation and washed with sterile 0.8% gelatin. Bacterial OD was adjusted to 5 X 
10 8 cfu/ml in 0.8% gelatin for detach leaf infection assay. Leaves of three week old 
transgenic lines and Col-0 were surface sterilized with 0.04% NaOCl and washed with 
autoclaved water several times. Finally, two spot of 10 ul bacterial inoculum was 
applied on each leaf and incubated at 28 C for 48 h to establish the infection and 
bacterial cfu was then determined. For inplanta spray assay, bacterial inoculum was 



adjusted to 5 X 10 5 cfu/ml in 0.8% gelatin and sprayed on whole plants equally. 
Infected plants were kept in growth chamber under 16 h light- 8 h dark cycle at 22°C 
for 7 days and then bacterial titres were measured. 

Sub-cellular localization. For sub-cellular localization, complete coding sequences of 
OsMATEl and OsMATE2 were fused upstream, in frame with, to sGFP 38 under 
control of CaMV 35S promoter in 326 sGFP vector 56 . These constructs were 
transiently expressed in protoplasts isolated from mesophyll cells of four weeks old 
Arabidopsis leaves and 15 u.g of purified DNA was used in polyethylene glycol 
mediated transformation 59 . After overnight incubation transient expression of GFP 
fusion protein was examined. Images were captured under a confocal microscope 
(LSM510, Zeiss, USA). 



Flavonoids estimation. For flavonoids quantification, individual flavonols in the 
rosette leaves were determined by methods of Misra et al. (2010). Briefly, plant 
material (500 mg) was grinded finely in liquid N 2 and extracted in 5 ml 
methanol : water (80 : 20) with constant agitation for 2 h and supernatant was 
collected after centrifugation. Pellet was re-extracted in 5 ml methanol : water for 
overnight at room temperature with constant agitation and supernatant was filtered. 
Both the filtrates were pooled, concentrated in rotavapour and again dissolved in 
1.0 ml of 80% methanol. Separation for qualitative and quantitative analysis of 
flavonols was performed by HPLC-PDA with a Shimadzu (Kyoto, Japan) LC-10 
system comprising an LC-10 AT dual pump system, an SPD-M20A PDA detector 
(254 nm) and rhodyne injection valve furnished with a 20 ul sample loop. 



Statistical analysis. Each experiment was carried out under a completely randomized 
design with at least three replications. Phenotypic observations were measured upto 
2-3 generations for consistency. Data are expressed as mean ± SEM and subjected to 
one-way ANOVA followed by post hoc Newman-Keuls for multiple comparison test 
of significance using GraphPad Prism 3.02. 
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